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We consider optimization of a rubidium atom clock that uses magnetically trapped Bose condensed 
f ^ ■ atoms in a highly elongated trap, and determine the optimal conditions for minimum Allan variance 

^ ' of the clock using microwave Ramsey fringe spectroscopy. Elimination of magnetic field shifts 

, and collisional shifts are considered. The effects of spin-dipolar relaxation are addressed in the 

optimization of the clock. We find that for the interstate interaction strength equal to or larger 
than the intrastate interaction strengths, a modulational instability results in phase separation and 
symmetry breaking of the two-component condensate composed of the ground and excited hyperfine 
clock levels, and this mechanism limits the clock accuracy. 
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PACS numbers: 95.55.Sh, 03.75.Kk, 03.75.Mn 



I. INTRODUCTION 



The most accurate atomic clocks in operation today are based either on trapped single ions or on atomic beams. 
The advantage of ion clocks is that a single ion can be trapped very tightly by static electric fields such that optical 
transitions do not cause significant heating or an escape of the ion. These optical transitions are decoupled from the 
trapping potential such that both ground and excited atomic states feel the same potential. On the other hand, one 
advantage of atomic beam clocks is the large number of atoms in a beam, such that the quantum projection noise 
can be decreased by 2-3 orders of magnitude with respect to the single ion clock. A clock based on a thermal atomic 
beam suffers from the velocity distribution that limits the transition bandwidth. Using a Bose-Einstein condensate 
fH ', (BEC) would significantly ameliorate this problem. A new type of atomic clock based on neutral atoms trapped in a 
O deep "magic-wavelength" optical lattice (magic because the transition does not have an optical light shift because the 
l— ~~ 1 , difference between the ac polarizabilities vanishes at the wavelength of the optical lattice) has recently been suggested 
^vq ' QLQ- This kind of clock can be operated on an optical transition, rather than a microwave transition, and promises 

| to be most accurate, but clocks of this type have not yet been fully characterized, 
■^j- ■ Here we consider whether a good atomic clock can be based on a more common type of trapped ultracold atom 
£SJ ' configuration, i.e., on a BEC in a magnetic trap. Trapped BECs can have many atoms, which gives them the large 
. number advantage mentioned above. This kind of clock can be much more accurate than a thermal cloud clock 
because the Doppler effect in a thermal clould can severely limit clock performance (see Sec. IIII B1 below). This effect 
is negligible for a BEC. In addition a BEC cloud has a well defined energy determined by the chemical potential 
that is uniform over the BEC and this helps in lowering the variance of the clock frequency. This kind of clock 
might be miniaturized, as microtraps for atomic BECs can be created above a fabricated chip. As has now been fully 
demonstrated, magnetic microchip traps can be used to manipulate neutral atoms on the micrometer scale Q. A 
high density, coherent atom source can be created via Bose-Einstein condensation on an atom chip Q , and "atomic 
conveyor belts", waveguides, and beam splitters can be implemented on atom chips Q. It is therefore intriguing to 
(-h ' entertain the possibility of creating an atomic clock on an atom chip . Hence, it is important to study theoretically 
O ,, and experimentally the potential of this kind of clock. One experiment of this kind, using Ramsey spectroscopy [31 , 
has already been carried out Q , and another experiment has been performed in a macroscopic magnetic trap @ , but 
using the same spectroscopic method. More recent mesoscopic atom clocks using coherent population trapping have 
been reported @ , and have the benefit of allowing compact optical light sources. 

Specifically, we consider a BEC in a magnetic trap and investigate a clock based on the Ramsey separated field 
spectroscopy method Q in a highly elongated trap. The quasi-lD geometry of an elongated trap has the advantage 
of further reducing the inelastic ultra-cold collisions as shown in Ref. 0]. As in Refs. @ and Q, we consider a two- 
photon microwave transition between two 87 Rb hyperfine states with an atomic frequency vq ~ 6.8 GHz. We treat the 
dynamics of the clock in mean-field and consider the amplitude and phase of the order parameters for the ground and 
excited clock states of the system, solving the coupled set of ID Gross-Pitaevskii equations to analyze the microwave 
clock frequency shift due to collisional and magnetic field effects. We determine the clock frequency shift introduced 
by the external magnetic potential and the kinetic energy of the Bose condensed gas, both of which are influenced by 
the difference in the size of the two atomic wave packets. The clock is designed to run with 87 Rb atoms in a magnetic 
field regime where the two hyperfine levels correlating with 5 1 Si/2\f = l,m/ = — f) and 5 1 5'i/2 1 / = 2,mj = f) 
experience the same first order Zeeman shift 0, [§|, [Io| • The collisional frequency shift from the resonance frequency 
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vn can be calculated d la Ref. [8j. As we shall see, the collisional shift can be cancelled by using the Zeeman shift 
Q and by optimizing the population difference in the ground and excited states [ll[. The latter is possible for an 
interstate interaction strength larger or smaller than both the ground and excited intrastate interaction strength. The 
clock run-time is limited by atom loss due to collisional spin dipolar collisional relaxation of the excited state 0]. 
For a 87 Rb condensate at high density (12.6 x 10 13 cm -3 ), the collisional dipolar loss in the excited state has been 
experimentally measured @. Atom loss led to a total density to drop by 3% in 20 ms. Nevertheless it is important 
to have a reasonable atomic number density n to compensate the effects of quantum fluctuations; the uncertainty, as 
quantified by the Allan standard deviation a scales as n, -1 / 2 . So an optimization of the density of atoms is necessary 
to reduce the quantum fluctuations and the collisional dipolar relaxation in order to increase the clock time. The first 
experiment with this type of clock using a trapped thermal cloud of 87 Rb atoms containing about 1.5 x 10 4 atoms 
with a density of atoms less than 5 x 10 12 cm~ 3 , yielded an Allan standard deviation 

<7(r) = 1.7 X lO" 11 ^ 1 / 2 Hz" 1 / 2 , (1) 

where r is the averaging time @. For the trap parameters used here, i.e., with radial frequency ui r j2-K = 120 Hz and 
axial frequency uiz/^tt = 0.5 Hz, the Allan standard deviation is of order of magnitude 10~ 12 y/T c /r, where T c is the 
cycle period. This trap confines the geometry to quasi-lD and has the advantage of further reducing the collisional 
dipolar relaxation Q. 

Some additional crucial limitations might make a magnetically trapped BEC unsuitable. A significant limitation 
is the collisional interaction between the atoms. Particularly problematic is the difference in the s-wave scattering 
length between atoms occupying different hyperfine levels which affects the collisional shift of the clock frequency. The 
collisional shifts of rubidium atoms are relatively small compared to cesium atoms [HI, HH, for example, but they can 
still be significant if many atoms are tightly trapped together. In order to minimize collisional shifts, we shall employ 
a method of overcoming collisional shifts by adjusting the ground to excited state ratio during the Ramsey fringe 
spectroscopy. Moreover, the run-time of the clock is also limited by the dynamics of the atomic cloud that can result 
in phase separation of the two spin components [Til , [ill . [r| [l7l . Ha ] . We find that a modulational instability results 
in the dynamics and the evolution depends on both the density of atoms and the balance between the interstate and 
intrastate interaction strengths. This gives rise to phase separation and symmetry breaking of the two-component 
condensate for the ground and excited clock levels that occurs after the first ir/2 Ramsey pulse that puts the atoms 
in a superposition of the ground and excited state. The modulation instability limits the clock accuracy. It therefore 
appears that magnetically trapped BEC clocks on an atom chip cannot promise to be the most accurate type of clock. 
The most significant limitation to the clock stability arises from the dynamics of the atomic cloud that creates a phase 
separation of the two wave packets for the ground and excited state. The time dependence of the phase separation 
depends on the density of atoms and on the interstate interaction strength; the smaller the density and/or the smaller 
the interstate interaction strength, the longer the phase separation time. Hence, a very weak axial trapping frequency 
(e.g., uj z /2tt < 0.5 Hz) resulting in a lower density of the atoms, allows an increased interrogation time and/or a 
greater total number of atoms, and therefore a further increase the stability of the clock beyond 10 _12 y / T c /r. 

The paper is organized as follows. The model of the clock based on Ramsey spectroscopy is described in Sec. UT1 
Section [ill Al briefly presents the numerical approach we use to analyze the clock. In Sec. lIIIBl we depict the quasi-lD 
dynamics in a trap that is very tight in two directions, and describe why the spin-relaxation collision mechanism, as 
well as other inelastic scattering processes, is suppressed in a ID geometry. Sec. IIII Cl describes the results obtained by 
numerically solving the coupled Gross-Pitaevskii equations for the order parameters of the ground and excited clock 
states. In Sec. IIVI we discuss two ways to improve and optimize the stability and accuracy of the clock by cancelling 
the collisional shift. Section IVl concludes the paper. 



II. MICROWAVE BEC MAGNETIC CLOCK USING RAMSEY FRINGES 



We consider an atomic BEC trapped in an external magnetic potential. The spatial variation is harmonic about 
the trap minimum. The atoms are initially in the ground electronic state, labeled |1), and a radio frequency field 
can transfer atoms into an excited state labeled |2). More specifically the levels \ f, m >= |2,1 > and |1, — 1 > arc 
used, and the transition involves a combination of a microwave pulse at 6.8 GHz to transfer the atoms from |2, 1 > 
to 1 1,0 > and then another RF pulse to transfer them from |1,0 > to |1, — 1 >. The |1, — 1 > state is trapped with 
the same potential as |2, 1 > if the magnetic field at the trap bottom is around 3.23 G [lCl ]. 

The clock described here uses the Ramsey separated field method 0- The atomic cloud interacts with two short 
microwave pulses separated by a time T; each pulse has pulse area close to 7r/2. A spatial inhomogeneity of the 
atomic energy levels is due to the spatially dependent Zeeman energy due to the magnetic field varying with position. 
Clearly, this can adversely affect the clock frequency. This effect is minimized by using a pair of energy levels which 
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experience the same trapping potential at a particular magnetic field strength. Refs. [5[ and [8j showed that at a 
magnetic field of ~ 3.23 G, the |1) = \ f = l,m/ = —1) and |2) = |/ = 2,m/ = 1) hyperfine levels of the 5Si/ 2 ground 
state of 87 Rb experience the same first order Zeeman shift such that the differential shift of the two levels across 
the cloud was ~ 1 Hz. The collisional shift also contributes to the spatial inhomogeneity of the atomic transition 
energy level across the cloud since the density of the cloud varies with position. However, as noted in Ref. [|[, it may 
be possible to cancel the Zeeman shift with the collisional shift. The stability and accuracy of the clock arc further 
improved and optimized by cancelling the collisional shift (as we shall see in Sec. IIVI) . 

The initial condensate starts in the ground state |1), and after the first ir/2 pulse, which we model by Bloch sphere 
dynamics assuming that the pulse duration r p is extremely fast compared to other time-scales, we solve a set of 
coupled Gross-Pitacvskii equations to describe the dynamics of the two-component (|1) and |2)) wave packets. After 
a time T a second short 7r/2 pulse is applied. For an intense short near-resonant pulse, the solutions of the optical 
Bloch equations for a two level atom gives the following unitary transformation operator: 

/ cos(ft g i/2) - i^- sin(fl g t/2) sin(fi a i/2) 
U W = 1 -i^sin(SV/2) cos(Q g i/2) 9 + sin(Q g t/2) 

where fl is the Rabi frequency, is the detuning and fl g = -\/|£!| 2 + Ais 2 is the generalized Rabi frequency. 
For example, if an atom is initially in state |1) and interacts with an on- resonant 7r/2 pulse, it evolves to the state 
(|1) + i\2))/y/2. This transformation can be used to describe the effects of both the first and second Ramsey pulses. 




III. MEAN-FIELD ANALYSIS OF CLOCK DYNAMICS 



The performance of the clock is affected by the dynamics of the two-component BEC after the first ir/2 Ramsey 
pulse. We shall see below that, because of the crossed interaction energy of the two spin components created after the 
first 7r/2 Ramsey pulse, the system becomes unstable, and the components eventually undergo a local phase separation 
that leads to symmetry breaking. The phase separation of the spin components limits the Ramsey interrogation time 
and hence the stability of the clock. 

In this section, we first describe the numerical methods used to investigate the clock dynamics in mean-field. Many- 
body effects can also be included as formulated in Ref. [2Cj, but we shall not do so here. Then we discuss the advantage 
of operating the clock in a highly elongated trap configuration. We present numerical results for this configuration 
and analyze them. 



A. Numerical method 



We investigate the clock dynamics in mean-field. The initial zero temperature condensate wave function (order 
parameter) is obtained by numerically determining the lowest eigenstate ip( r ) of the time-dependent Gross-Pitaevskii 
equation for particles of mass m, confined in an external potential V cx t(j) and a mean-field interaction energy due to 
contact two-body interactions with coupling strength gn = 4Trh 2 an/m where an is the s-wave scattering length for 
atoms in the ground state. This is accomplished with an imaginary time split-step Fourier transform method. The 
effect of the first pulse that couples the two atomic spin states is modeled using a unitary transformation on the zero 
temperature ground state wave function and gives two wave functions representing the ground state and the excited 
state atoms: 

^(r,0) = A^(r) (3) 

where £ = 1,2 correspond to ground and excited state labels respectively, and Aj is the complex amplitude of state i 
obtained using Eq. ([2|). We take the normalization of the initial condensate wave function such that / \i>(r, 0) | dr = N, 
where N is the total number of atoms, and the amplitudes Aj are determined by the Bloch sphere dynamics for the 
two levels in the presence of the microwave field inducing the transition . The amplitudes Ai satisfy < | Aj | < 1 
and MM = 1- The two component condensates evolve according to 

ift*Wp2 = - + V- cxt (r) + (-l)^il^L + ^ = i, 2 ^|^(r, 0| 2 ^)^(r,*) (4) 

where the atomic resonance transition frequency is denoted asw 2 i- The interaction strength g^- = Anf^a^ a,ij /to, with 
i,j = 1,2, is defined in terms of the s-wave scattering length for particles in states i and j, aij, and the two-particle 
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correlation parameter for particles in states i and j at zero separation between particles, oq) pi. I2H I22I, I23I, l24l, [25| . 
The latter quantity is often denoted as . The values of the two-particle correlation parameter is such that 

(2) (2) (2) 

< < 2 and 1 < ■ ^ ■ < 2 for Bosons. For a condensate, = 1. For the inter-state (two-component) 

(2) . . (2) (2) 

two-particle correlation at zero separation, a 12 in a condensate we considered two values, a 12 = 1 an d a 12 = 2. 
We have learned recently that was measured to be nearly unity for 87 Rb [2(|, but to we shall present results of 
calculations for both values. 

We propagate the two BEC components for a time T between the Ramsey pulses by solving the coupled time- 
dependent Gross-Pitaevskii equations using the split step Fourier transform method. During the propagation for a 
time T, the phase and amplitude of the component wave functions evolve with time. After the time T, we apply again 
the unitary transformation operator corresponding to a tt/2 pulse. We then integrate the component wave functions 
over space to obtain the probabilities for finding the atoms in the two states. 



B. Dynamics of the clock in a highly elongated trap 

The effects of the mean-field collisional dipolar relaxation are important to investigate in order to optimize the 
accuracy and the stability of the clock. As we shall see, phase separation of the two spin components due to mean- 
field dynamics, and loss of excited state atoms due to the collisional dipolar relaxation between atoms in the excited 
state arc two factors that can significantly reduce the performance of the clock. The collisional dipolar loss of the 
excited state for a 87 Rb condensate at high density (12.6 x 10 13 cm -3 ) has been experimentally measured Q. Atom 
loss caused the total density to drop by 3% in 20 ms. In order to reduce this effect, we propose to run the clock in a 
highly elongated trap; in a quasi-one-dimensional condensate the spin-dipole relaxation collisional loss can be made 
much smaller than in 3D Q. The inelastic rate coefficient for going from the incident channel f3 to a final inelastic 
channel /3 for strong confinement K con{ p> p is related to the 3D inelastic rate K free p/p by 0], 



Here po is the relative collision momentum and E = j>Q/2/i r is the relative collision energy in the incident channel, 
Ppi = yj2fi r (Eo + Epi p) is the relative momentum in the final inelastic /3 channel, where Epi p is the asymptotic 
energy difference of the two channels, fi r = m/2 is the reduced mass, and the parameter Up>p is the coupling between 

the channels (3 and j3. This quasi-lD rate is reduced by a factor of jW^bi from the 3D rate, and is small for 
ultra-cold collision energies. For high anisotropy of the trap, the ID interaction strength is given by 

(6) 

J ma i 



where a± = ^Jhjmui r is the radial harmonic oscillator length and Lo r /2n the radial trap frequency [27] . For sufficiently 
large u) r , radial profile is harmonic oscillator like, and the motion of the atoms are frozen in the radial direction. We 
assume that the magnetic field is such that the first order Zeeman shift is the same for the two atomic internal states. 
Now, if the frequency uj z /2it is also made small, so the ID density is small, the nonlinear interaction term gjfn can 
be made very small. 

Using a magnetically trapped thermal cloud can significantly reduce the mean-field collisional shift compared to a 
BEC, however, a thermal cloud has a Dopplcr width that increases the bandwidth of the clock transition and can 
therefore limit the clock accuracy. Indeed the clock bandwidth at half maximum is in our case inversely proportional 
to the interrogation time T; the larger the interrogation time, the smaller the transition bandwith of the clock, and 
the Doppler effect on the bandwidth becomes more significant. With our interrogation time ( 0.5 s), the Dopplcr 
width must be much smaller than 1 Hz to be negligible. Hence, the temperature of the thermal cloud must be well 
below 2.5 /xK so as to be competitive with a BEC. For a 500 nK 87 Rb thermal cloud (i.e., above condensation) the 
Dopplcr width is 0.4 Hz which is smaller than the bandwith of the clock transition. Moreover, in a BEC the energy 
is determined by the chemical potential that is uniform across the BEC and very well specified; this helps in lowering 
the variance of the clock frequency. The collisional frequency shift calculated with 10 4 atoms in an highly elongated 
BEC corresponding to a density of 3 x 10 14 atoms/cm 3 for our trap geometry, is about 3 Hz, but the uncertainty in 
this collisional shift (due to the uncertainty in the number of atoms and the uncertainty in the value of the scattering 



5 




-5 5 



10 15 



-15 -10 -5 5 

z[aj 



10 15 



FIG. 1: (color online) Phase of condensate ground state (solid curve) and excited state (dashed curve) components for the 
interstate two-particle correlation parameter at zero separation ot^ = 2, after the first 7r/2 pulse as a function of position, z, 
in the magnetic trap at times (a) t = ms, (b) t = 6.3 ms, (c) t = 180 ms and (d) t = 500 ms. a z = \Jh/muj z is the axial 
harmonic oscillator length. 



length) is very much smaller than 3 Hz, and hence the collisional shift can be largely compensated for as far as the 
clock frequency is concerned. Although the collisional shift for a 500 nK thermal cloud is at least 100 times smaller 
than that for the BEC, the Doppler width of the thermal cloud is 0.4 Hz, and this is presumably much larger than 
the uncertainty in the collisional shift of the BEC. 



C. Numerical Results 



Our numerical calculations have been carried out with an axial (y z ) and radial (v r ) trap frequency of 0.5 Hz 
and 120 Hz respectively, so the anisotropy ratio A = v r jv z = 240. The three scattering lengths are taken to be 
an = 100.44 ao, «22 = 95.47 ao and a\i = 98.09 ao Q, where ao is the Bohr radius. In order to better understand the 

(2) (2) 

effect of the modulational instability on the clock, we carry out the calculation with a 11 22 = 1 and a 12 = 2. The 
optimization of the number of atoms for the frequencies given above that gives the best Allan deviation is of the order 
of magnitude 10 atoms. The two pulses used for the Ramsey separated field method are taken to be 7r/2 pulses, i.e., 
\A,\ 2 = \A 2 \ 2 = 1/2. 

Fig. [1] shows the phase 9i(z,t) of the i th condensate wave function, %pi(z,t) = |^>i(z,i)| exp[i6i(z, i)], as a function 

(2) 

of position, z, in the magnetic trap for a 12 = 2. Immediately after the first 7r/2 pulse, the phase of the two spin 
components is spatially uniform and their difference is 7r/2 (Fig.[T](a)) as is easily understood from the transformation 
in Eq. (J3J). Following the 7r/2 pulse, mean-field effects begin to create a spatially varying phase across the two 
condensate wave packets (Fig. [T](b) and (c)) (28|. Beyond t = 0.24 s, the spatially dependent variations in the phase 
appears completely chaotic; the mean-field treatment has not only reached the point of numerically limited accuracy 
but has actually lost its regime of validity. 

Fig- HI shows the evolution of the position dependent density of the two atomic states for a sequence of inter-pulse 
times T = s (Fig. HJa)), T = 0.18 s (Fig. HJb)), T = 0.22 s (Fig. [^c)), T = 0.26 s (Fig. [DJd)), T = 0.5 s (Fig. [!»), 

and T = 1 s (Fig. [2jf ) ) , and a 12 = 2. The density of the two components are almost identical at T = s and 
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FIG. 2: (color online) Condensate density of ground state (solid curve) and excited state (dashed curve) components for the 
interstate two-particle correlation parameter at zero separation = 2, after the first 7r/2 pulse for a sequence of inter-pulse 
times (a) T = s, (b) T = 0.18 s, (c) T = 0.22 s, (d) T = 0.26 s, (e) T = 0.5 s and (f) T = 1 s. At T = s the 
two spin components have the same amplitude since the microwave pulse is very short but dephased of 7r/2 for a n/2 pulse. 
a z = ^Jh/muiz is the axial harmonic oscillator length. 



T = 0.18 s and are smoothly varying with position. The density profiles begins becoming irregular at the center 
of the condensates at about T = 0.20 s. At later times a spiked structure whose amplitude increases with time 
develops and some local phase separation occurs (Fig. E3c)-(f)) due to the repulsive interaction between the wave 
packet components. A similar spiked structure of the density as a function of position has been obtained numerically 
in the regime of strong excitation of the BEC loaded in a ID optical lattice plus an asymmetric external magnetic trap 
by instantaneously giving a large displacement to the initial position of the center of the magnetic trap in Ref . (29| . 

The dynamics is different for a 12 — 1- Fig. [3] shows the evolution of the position dependent density of the two 
atomic states for a sequence of inter-pulse times T = s (Fig. IHa)), T = 0.36 s (Fig. H^b)), T = 1 s (Fig.(3Jc)), T = 2 
s (Fig- EUd)), T = 50 s (Fig. E^e)), and T = 220 s (Fig. [3Jf)). The phase separation appears later, around 0.4 s, and 
evolves more slowly than the case of strong interstate interaction strength (a 12 = 2). Since the interaction strength 
of the ground state is larger than in the excited state, the ground state density protrudes beyond the excited state 
density (T = 1 s), and at later times the phase separation evolves in a complex way under the effect of the instability 
(T = 2 s and 50 s) and eventually a symmetry breaking occurs (T = 220 s) . 

A rough estimate of the time scale at which the system becomes sensitive to the phase-separation instability, T ps , 
can be obtained by using the expression derived by Timmermans in Ref. (l7j for a homogeneous system: 

r ps = 27r/|n_ lfc/ 1 = 2Trh/m\c. | 2 , (7) 

where 



ft 2 _ ik = c 2 _k 2 + (hk 2 /2m) 2 (8) 
is the dispersion of the double condensate excitation, and the parameter 

h 2 



c 2 



a^aunx + a 2 2 ^a 2 2n 2 - \J (a^ anni) 2 + {a^a 22 n 2 ) 2 + (2(a^ ) a 12 ) 2 - a 2 2 2 an^^n^ J (9) 
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FIG. 3: (color online) Condensate density of ground state (solid curve) and excited state (dashed curve) components for the 
interstate two-particle correlation parameter at zero separation = 1, after the first n/2 pulse for a sequence of inter-pulse 
times (a) T = s, (b) T = 0.36 s, (c) T = 1 s, (d) T = 2 s, (e) T = 50 s and (f) T = 220 s. At T = s the two spin components 
have the same amplitude since the microwave pulse is very short but dephased of n/2 for a n/2 pulse. 



is the phonon-like sound velocity at low momenta. For instability, c_ < and f2_ fc < 0. The fastest growing mode 

has wave number kf = v2m\c-\/fr, and grows with an initial rate of m\c-\ 2 /h. Eqs. ©-© show that the local phase 
separation of the two condensate components and the symmetry breaking are due to the cross interaction terms under 
the condition c_ < which occurs when (a 12 a i2) 2 > a u "22 111022/2. Eq. © shows that the time at which the 
symmetry breaking starts depends on the atomic density; the higher the density (corresponding to larger parameter 
|c_|), the smaller the symmetry breaking appearance time. We obtain a value of r ps = 35 ms which is smaller than 
the time from which the phase separation starts in Fig. [51 This discrepancy is probably due to the space-dependent 
atomic density that locally changes the value of r ps , increasing from the center of the clouds where the density is 
higher to the edge where the density is smaller. As the instability depends on the density of atoms, the local spiked 

structures start to develop at the center of the trap and then spread throughout the clouds. 

(2) 

Note that the case of a 12 = 1 with the same numerical parameters values yields T ps « 7 s and the numerical 
calculation yields a value beyond 120 s for the appearance of the symmetry breaking (Fig. [6]). In this case, the phase 
separation appears much earlier, w 0.4 s before the symmetry breaking takes place. Indeed as the intrastate and 
interstate interaction strengths arc almost identical, the BEC clouds are not strongly perturbed after the first ir/2 

pulse. Hence, the two BEC components are barely unstable and phase separation and the symmetry breaking does 

(2) 

not occur until much later than in the case of a\ 2 ' = 2. This clearly increases the stability of the clock. The smaller 
the difference between the interstate and intrastate interaction strength, the longer the interrogation time and the 
higher the stability of the clock. 

The spiked structure and phase separation in position space is correlated with a derealization in momentum space 
created by the strong excitation of the two BEC components as the dynamics proceeds. Indeed the calculated density 
of the atoms in momentum space is completely dclocalized beyond 0.5 s, as shown in Fig. 2] The increased width of 
the momentum distribution observed from t = s to t = 0.22 s is due to the fact that after the first 7r/2 pulse the 
interaction energy converts to the kinetic energy. To monitor the symmetry breaking shown in Fig. [TJd) , we calculate 

the mean value of the axial momentum of each spin component as a function of time, (pi z (i)) = J Pi z {t)\ipi(p z , t)\ 2 dp z 
where ipi(jp z ,t) is the Fourier transform of the condensate wave function tpi(z,t). Fig. [5] shows the evolution of the 
mean value of momentum of the two condensates as a function of time. The mean value of the momentum for each 
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FIG. 4: (color online) Condensate momentum density of ground (solid curve) and excited (dashed curve) state component for 
the interstate two-particle correlation parameter at zero separation ct^ = 2, after the first 7r/2 pulse for a sequence of times 
(a) T = s, (b) T = 0.18 s, (c) T = 0.22 s, (d) T = 0.26 s, (e) T = 0.5 s and (f) T = 1 s. 



spin component is zero from to 0.24 s, and then the mean momentum of each component starts oscillating in time, 
but the total momentum is conserved. We note that the symmetry breaking appears after the density starts showing 
an irregular profile of small amplitude at its center. Comparison of Fig. [5] and Fig. [6] shows that phase separation 
occurs at much later times. 

We now analyze the effect of the mean-field dynamics of the system on the performance of the clock. Fig. [7] shows 
the calculated probability P e for finding atoms in the excited state immediately after the second 7r/2 Ramsey pulse as 
a function of the detuning Ais of the microwave frequency from the atomic transition frequency times the interrogation 
time, /S.vT. The curves in Fig. Ba) have been calculated for three values of the interrogation time, T = 0.18 s (solid 
curve), 0.5 s (dashed curve), and 1 s (dotted curve). The fringe contrast decreases as the interrogation time increases. 
Fig.[7[b) plots the variance of the excited population NP e (l — P e ) [I3,|3l]]. The variance is large at time T = 0.5 and 
1.0 s where the condensate profiles are spiked and asymmetrical, and this results in poor stability of the clock. The 
interrogation time of T = 0.18 s where the condensates show a smooth and symmetrical profile lead to a frequency 
stability of 2.6 x 10~ 12 yjT c /r where T c is the cycle period and r the averaging time. It is important to note that 
the stability can be improved by further decreasing the axial and radial frequencies, keeping a high anisotropy ratio 
so that the quasi-lD regime remains. The improvement will depend on how low the axial frequency can be made 
without causing fluctuations of the trapping magnetic field. The goal is to further lower the density so as to increase 
the time at which the phase separation of the two spin components after the first microwave pulse, and hence increase 
the interrogation time. The total number of atoms can also be independently optimized. Thus, it will hopefully be 
possible to reach a stability beyond 10 _12 y / T c /r. 

We can express the Allan standard deviation for a Ramsey fringe experiment as a = — ^Vn VTc/t, as a function 
of the axial trapping frequency. The interrogation time is fixed by the density of the atoms; the higher the density, 
the smaller the interrogation time. By fixing the interrogation time Tfi x to correspond to a given density of atoms 
(n = N/lg, where l z = (3^f iV/2mw 2 ) 1/3 it is possible to vary the number of atoms and the axial frequency. By 
equating the density calculated for different number of atoms and trap frequency, we derive the number of the atoms 
as function of the axial trapping frequency: 



(10) 
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FIG. 5: (color online) Evolution of the mean value of the momentum of ground (solid curve) and excited (dashed curve) state 
component after the first 7r/2 pulse as a function of time for ot^ = 2. 



where A^ x and uj z ,&x are respectively the given number of atoms and axial trapping frequency that fix the density 
and the interrogation time Tg x . Then the Allan standard deviation is: 

1 / U) t \l/2 



and varies as the square root of the axial trapping frequency. In our case, using iVg x = 10 4 atoms and w Zj fi x = 0.5 

(2) 

Hz, we find Tg x = 180 ms for a 12 — 2, and the new values of the interrogation time obtained by varying the axial 
trapping frequency arc within about 10% of Tfi x . 

The projection noise does not have a simple expression as a function of the radial frequency because the instability 
of the BEC depends both on the density and the geometry of the cloud. Furthermore, the variation of the radial 
frequency deeply affects the geometry of the cloud and changes the interaction strength g\f . 

Normally we should add to the projection noise, the frequency noise due to the fluctuation of the density introduced 
by the fluctuation of the total number of atoms. However, as we shall sec in the next section, there is a simple method 
to cancel the collisional frequency shift and then the noise due to density fluctuations will cancel too. 



IV. IMPROVEMENT OF THE CLOCK: CANCELLATION OF THE COLLISIONAL SHIFT 

In this section we analyze the possibility of improving the clock by cancelling the collisional shift. One of the 
advantages resulting from the cancellation of the collisional shift is that the clock becomes insensitive to the variation 
of the total number of atoms. Indeed, even if the density is low, nevertheless the variation of the total number of 
atoms at each cycle period creates a variation of the density of the atoms and this introduces noise that limits the 
stability of the clock. As the collisional frequency shift depends on the density of atoms, the idea is to cancel the 
collisional shift so that such a variation of the density of atoms does not affect the stability of the clock. Note that the 
variation of the number of atoms we are referring to is not due to quantum fluctuations but rather due to experimental 
fluctuations in the number of atoms. Here we present two different proposals to overcome this problem. 

A small density of atoms in the excited state minimizes the collisional dipolar relaxation loss. But the clock 
sensitivity to quantum fluctuations increases if the population of atoms in the excited state is too small. So an 
optimization of the density of atoms is necessary to obtain both good signal-to-noise ratio and a long clock time, yet 
having small collisional dipolar relaxation. We can use the Zeeman shift to compensate the collisional shift. This 
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FIG. 6: (color online) Evolution of the mean value of the momentum of ground (solid curve) and excited (dashed curve) state 
component after the first n/2 pulse as a function of time for ct^ = 1- 



method can be applied both to a normal cold atomic cloud and a BEC, and has been proposed in Ref. [8[. However, 
this method requires adjusting two parameters and seems not to be simple to implement. We will not use it in this 
paper. Instead, we study the cancellation of the collisional shift terms by playing them off against each other as 
discussed by Gibble and Verhaar in Ref. [ll[ for a thermal cesium atom clock. The expression for the collisional 
frequency shift is |, ED, SI M, [H 

A ^ I ( 2 ) i ( 2 ) ( 2 ) (2) \ Mn\ 

Aj/i n t = 2"( a i2 a i2"4 + «22 "22^2 - ol\\ a\\n\ - a{ 2 a 12 n 2 ) ■ (12) 

TTITT CL j 

where n\ and n 2 are the density of atoms per unit length in the ground and excited state respectively By equating 
Ai/int to zero, wc obtain a simple relation between the density of the atoms in the two states: 

(2) (2) 

n 2 _ a\{ai2 - a. x {a\\ . . 

m - (2) (2) ' [ > 

"i a\ 2 'a\2 - a 22 'a 2 2 

The above relation is satisfied only if a^ 2 > otyy an/012 and ct[ 2 > a^ 22 a 2 2/ai2; or a^ 2 < ctii/ai2 and < 

f2) f2l f2) ^ " 

a 2 2 "22/012 • When = = 1, and 0,12 is only slightly different from an and 022, such that 0,2 < an < an, 

(2) (2) (2) 

the cancellation of the collisional shift is possible only for a 12 > 1.02 or a 12 < 0.97. For = 1 the ratio n 2 jn\ in 
Eq. (fT"3|) is negative, and therefore cancellation of the collisional shift is not possible. For the scattering lengths of 

87 Rb and a^ 2 = 2, the density ratio is n 2 = 0.95 n\. Such a ratio should maintain the sizes of the two components 
to be almost identical. For a ratio of densities close to unity, the Zeeman frequency shift terms due to the size of the 
atomic clouds will be small. However, a problem originates from |2)-|2) collisional dipolar relaxation which can be 
important at high density and this can limit the clock run-time. But in a quasi-lD system, the inelastic ultra-cold 
collisions are further reduced. To illustrate the improvement of collisional shift, Fig. [8] shows the calculated excited 
state population as a function of AvT for an inter-pulse time T = 0.18 s. The dashed curve is for the case of 7r/2 
pulses and the solid curve shows the case of the slightly different populations suggested above. In the latter case, the 
collisional shift is cancelled and this gives a better fringe contrast and a smaller frequency shift. The improvement 
of fringe amplitude is due to the fact that the difference in the phase of the condensates is space-independent as the 
collisional shift is cancelled. 
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The remaining frequency shift is now due to kinetic and Zecman shift. The frequency shift due to the kinetic term, 



Akinetic - , (14) 

is —10.997 Hz and the shift due to the Zeeman terms, 

(mu 2 z z 2 ) 2 - (mulz 2 )! 

^Zeeman — — T , (,-LOJ 

Ann 

is —0.266 Hz, giving a total frequency shift of Av = —11.263 Hz, in good agreement with the frequency shift 11.12 Hz 
determined from the interference pattern of the solid curve in Fig. [S] In Eqs. (|14j|15p , the symbol (. . denotes the 
expectation value calculated with the wave function ipi (t = T). We see that the shift introduced by the Zeeman terms 
is small because it is proportional to the difference of the square of the size of each condensate and this difference is 
small. The shift due to the kinetic terms is proportional to the difference of the square of the width of the spectral 
density of each condensate. As the size of the condensates and the width of their spectral density change in time, the 
shifts arc time dependent. 



V. CONCLUSION 



We modeled a microwave frequency atomic clock using a configuration of BEC atoms in a highly elongated magnetic 
trap. We showed that the stability of the clock for a trap radial frequency uj r /2-K = 120 Hz and axial frequency 
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AvT 

FIG. 8: (color online) Excited state population as a function of the detuning Ai/ of the microwave pulse from atomic transition 
times the interrogation time in the presence (dashed curve) or without collisional (solid curve) for an inter-pulse time T = 0.18 

(2) 

s, and the interstate two-particle correlation at zero separation a 12 = 2. 



lu z /2it = 0.5 Hz is 2.6 x 10 _12 y / T^/V and 1 x 10~ 12 \/T c /t for = 2 and afj = 1 respectively. The performance 
of the clock is related to the configuration of the trap and can be improved by running the clock with an even weaker 
axial trapping frequency. We found a dynamical instability that results in phase separation and limits the clock 
stability and accuracy if a long interrogation time is used depending on the density of atoms and on the ratio between 
the intrastate and interstate two-particle correlation parameter at zero separation. We considered optimization of 
the experimental parameters maximize the stability and accuracy of the clock. For a 87 Rb BEC, the collisional shift 
terms can be cancelled by playing them off against each other by adjusting the population in the ground and excited 
states with the Ramsey pulses so that the problem of density fluctuations does not strongly affect the stability and 
the accuracy of the clock. This can be done only if the two-particle correlation parameters satisfy the relations 

Oiy} > Oi\^a\\j a\2 and > c4^ a 22/ a i2, or < an/aiz and < a 22' ' a 22/ '«i2- The equality of these 
previous relations automatically cancel the collisional shift. 
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